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A SIMPLE INTEGRAL METHOD FOR THE CALCULATION OF 
REAL-GAS TURBULENT BOUNDARY LAYERS 
WITH VARIABLE EDGE ENTROPY 

By Charles B. Johnson and Lillian R. Boney 
Langley Research Center 

SUMMARY 

A simple integral method for the calculation of an equilibrium air turbulent bound- 
ary layer with variable entropy effects on a blunt axisymmetric body is presented. A 
modified Crocco enthalpy velocity relationship is used for the enthalpy profiles and an 
empirical correlation of the N-power law profile is used for the velocity profiles. The 
skin-friction coefficient expressions of Van Driest and Spalding and Chi, corrected for 
axisymmetric flow by a turbulent Mangier transformation are used in the solution of the 
momentum equation. Heat-transfer predictions are obtained by use of various modified 
forms of Reynolds analogy. 

Experimental heat- transfer data from wind-tunnel tests on a flat plate at Mach 6.0 
and on blunt 5° half-angle cones at Mach 8.0 and 11.5 are compared with the results 
of the method and good agreement is obtained. In addition, sample calculations are 
presented for a 5® half-angle cone at flight conditions of approximately 6.09 km/sec 
(20 000 ft/sec) at an altitude of 18.29 km (60 000 ft). The calculations show the effect 
of the variable entropy on the flow properties external to the boundary layer as well as 
on the boundary- layer integral parameters. 

INTRODUCTION 

In recent years the ability to predict the properties of compressible turbulent 
boundary layers has been developed to a significant degree as evidenced by the numerous 
sophisticated finite-difference methods of solution. A review of most of the compressible 
turbulent boundary-layer-calculation methods prior to 1967 is given in references 1 to 3 
and the significant developments in computation methods since 1967 are given in refer- 
ence 4. Integral methods have been used to compute the turbulent boundary layer for 
many years and these methods have some distinct advantages over the more sophisticated 
finite-difference methods. In some instances, in reference 5, for example, it was found 
there was no significant gain in accuracy for a particular finite-difference numerical 



solution compared with the three less rigorous, but much faster, integral methods. Fur- 
thermore, when real-gas and variable -entropy considerations are added to a turbulent- 
boundary-layer solution, the simplicity of an integral method becomes even more 
desirable . 

This paper presents a simple integral method for computing the compressible tur- 
bulent boundary layer with equilibrium air chemistry and variable-entropy edge effects. 
This paper describes in some detail the method and techniques used in the computer pro- 
gram to make the transitional and turbulent boundary- layer calculations. The method 
uses velocity and enthalpy profile correlations which are based directly on ground-test 
experimental data. The calculations are carried out in the physical plane, by using either 
the Van Driest n (ref. 6) or the Spalding- Chi (ref. 7) skin-friction relations and the 
Reynolds analogy to obtain heat transfer. The method is applied herein to a blunt axlsym- 
metric body; however, the method has shown good agreement with data for nozzle wall 
boundary layers over a wide range of conditions. (See ref. 8.) 

SYMBOLS 


Values are given in both SI and U.S. Customary Units. The measurements and cal- 
culations were made in U.S. Customary Units. 

A area 

ArjBr coefficients in equations (10) and (12) 

a speed of sound 

Cp average skin-friction coefficient based on conditions at edge of boundary layer 

Cf local skin-friction coefficient based on conditions at edge of boundary layer 

functions given by equations (20), (21), and (22) 

Fmt Mangier transformation factor 

g^ a dimensional constant 

H total enthalpy 

h static enthalpy 
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h heat-transfer coefficient 

M Mach number 

N exponent in velocity-profile relation 

Npr Prandtl number 

Nr recovery factor 

^St,e local Stanton number based on conditions at the edge of the boundary layer 

p pressure 

q normalized heat-transfer rate 

q heating rate 

R unit Reynolds number 

Raf Reynolds analogy factor 

Re,x local Reynolds number based on surface distance 

Rg^0 local Reynolds number based on momentum thickness 

dimensionless body radius 
r^ nose radius 

^s/^n dimensionless radius out to shock (entering stream tube) 

S entropy 

T temperature 

u velocity along x/r^ 

V velocity along y/r^ 
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X 


surface distance 


xy'rn normalized surface distance from stagnation point or sharp- cone apex 

XyQ distance from virtual origin of turbulent boundary layer 

y/r^ normalized coordinate normal to wall 

z correlation parameter defined by equation (8) 

Q!j^ exponent in equation (10) 

y ratio of specific heats 

6/rn dimensionless boundary -layer thickness 

6*/rjj dimensionless displacement thickness 

gyrjj dimensionless momentum thickness 

cone half -angle 
p density 

T shear stress 

Subscripts: 

aw adiabatic wall 

e edge of boundary layer 

i incompressible value 

max maximum 

R one of three regions in boundary layer, I, n, or HI 

s shock conditions 
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t stagnation conditions 

tr transition 

w wall 

«> free- stream value 


Primes denote evaluation at reference enthalpy condition of equations (13). 

DESCRIPTION OF THEORY 


An integral method for the calculation of the equilibrium air compressible turbulent 
boundary layer over blunt cones including the effects of variable entropy is presented. 

The present method has been developed from a revision of the integral method described 
in reference 9. The new method has also been used extensively for nozzle boundary- 
layer calculations over a wide range of conditions. Results are compared with data on 
nozzle walls at Mach numbers of 6, 8, and 19 in reference 8. 


Variable-Entropy Momentum Integral Equation 

The boundary- layer equations for the conservation of mass and momentum for 
application to a body of revolution at zero angle of attack are 


8pu ^ 
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where 


P = P(x) 


and 



Equations (1) and (2) are combined and integrated from the wall to the edge of the 
boundary layer. The resulting equations are cast in an integral form amenable for 
numerical integration as 
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Equation (3) is similar to the well-known momentum integral equation (eq. (42), Ch. 9 of 
ref. 10), with the exception of the additional term 
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For configurations with slight nose bluntness, a variable -entropy condition exists along 
the edge of the boundary layer. Therefore, the Bernoulli equation is not applicable and 
the previous term (eq. (6)) is added. 


Velocity Profile 

The velocity profiles used in the integral parameters (eqs. (4) and (5)) were calcu- 
lated from the N-power law relations: 


u 

ue 


/y/rn\^/^ 

Ivw 


(7) 


The value of N was calculated from a correlation of N which is a function of 

Rg a, Mg, T^/Tg, and shown in figure 1 . The limited amount of data for N 

’ ' ^n/ 

values used in the correlation were taken from turbulent velocity profiles on flat plates, 
cones, hollow cylinders, and shock tube walls taken from references 11 to 24 and corre- 
late as 

N = 6.0 log z - 7.0 

where 
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A small change in the value of N (that is, approximately ±10 percent) does not have a 
large effect on the calculated values of the integral parameters, and S/r^, pro- 

vided the value of N is greater than approximately 7. In general, the overall result of 
the boundary-layer calculation (particularly heat transfer) is insensitive to small changes 
in the value of N provided the velocity profile is fully turbulent (N > 7) . 


The parameter used in the z correlating parameter (eq. (8)) is similar 

^n/ ^n 

to the x/6 parameter used in the "N-power-law" correlating parameter of reference 25. 
The Xyo/^n used in equation (8) is the dimensionless distance from the virtual origin of 
turbulent flow, whereas the x used in reference 25 is the distance along the nozzle cen- 
ter line from the throat to the measuring station. The expression used to find the virtual 
origin distance on a cone is 



Cxp 

1.045 

2 


(9) 


where Cp^2 is the Spalding-Chi (ref. 7) average value of skin friction and is the 

local value of momentum thickness. The flat plate Cjj/2 from reference 7, used in 
equation (9), is correlated for use on a cone by multiplying the C^2 by 1.045. (See 
appendix.) The Xyg/r^ in equation (8) is always calculated from the Spalding-Chi (ref. 7) 
Cpy/2 regardless of what skin-friction law is used in equation (3), in order to be consis- 
tent with the values of Xyo/rn that were used to establish the correlation in figure 1. 

The exponents of the parameters in z (eq. (8)) were evaluated on a trial-and-error basis 
for over 300 combinations until the best combination of correlating parameters was found. 
A more detailed survey of the "N power-law" velocity profiles is given in reference 26. 


Density Profiles 


Boundary- layer density profiles were calculated for equilibrium air by use of the 
same real- gas thermodynamic tables that were used in the inviscid calculation. (See 
ref. 27.) The density at a point in the boundary layer was found by entering the equilib- 
rium air thermodynamic tables with the local static pressure and the calculated value of 
static enthalpy. The calculation of a density profile through the boundary layer is based 
on the assumption that the total enthalpy through the boundary layer varies as a function 
of u/Ug according to a modified form of the Crocco expression given by 


H-hw 
He - 


= Ap 



( 10 ) 


The total enthalpy at any point in the boundary layer is given by 
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H = h + 1^ (11) 

When equations (10) and (11) are combined, the resulting expression for static enthalpy is 


h = 


hw 




^e^ /u 

2 (uj 


( 12 ) 


where the coefficients Ar and Bfj and the exponent o!pj vary with three regions of 
calculation through the boundary layer (denoted by subscripts I, n, and HI) which are 
as follows: 


Region I (wall) (O ^ u/Ug < O.Ol).- The coefficients for the wall region as used in 
equation (12) are 

Ai = 0 I 


_ Haw ~ (^P^)w I 

I He - h-^ , .2/3 / 

(^Prj 

^1=1.0 


(13) 


where the coefficient Bj is derived from 



(14) 


combined with the energy equation, the Fourier conduction law, and the Colburn form of 
the Reynolds analogy. (See ref. 28.) The prime indicates that the Prandtl number is 
evaluated at Eckert's reference enthalpy which is 


h’ = i(hw + he) + 0 .22(haw - he) (15) 

where haw = he + Nr(He - he) and N^. = 0.89 were used. The local Npj. is found 
from table IV of reference 29 where the Prandtl number is tabulated as a function of tem- 
perature and pressures. The temperature used to determine the Prandtl number in ref- 
erence 29 is found from the real-gas thermodynamic tables as a function of h' (eq. (15)) 
and the local pressure pg. 

Region III (outer) (0.1 < u/ue = l.O).- The coefficients for the outer region (region in) 
as used in equation (10) are 
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Am = 0 

Bjn = 1.0 

ani = Constant at a given x station 


( 16 ) 




(Note that Q!]jj = 1.0 for sharp flat plate and sharp or slightly blunt cones and qjjjj = 2.0 
for axisymmetric nozzles. (See ref. 8.)) 


Region II (intermediate) (O.Ol ^ u/Ug g 0.1^.- Redons I and III are joined by an inter- 
mediate linear relationship which matches the region I (wall) equation at uy'Ug = 0.01 and 
the region HI (outer) at u/Ug = 0.10, where the coefficients used in equation (10) for the 
intermediate region are calculated from 

<^IIT 

Afi = Bjjj(O.l) - O.lBjj 


a 


Bii = 


Bm(O-i) 


m 


Bi(O.Ol) 


(0.09) 


(17) 


“11 = 


1.0 


J 


It should be noted that the boundaries of the three regions may be different from the val- 
ues assumed herein. (See ref. 8.) 


Mass Flow Into Boundary Layer 


The fact that the Bernoulli equation cannot be used in the momentum equation 
because of a change in entropy along the edge of the boundary layer has been previously 
pointed out in the derivation of equation (3) . There is a change in entropy between the two 
closely spaced streamlines passing through the shock and this change in entropy at the 
shock is transmitted to the edge of the boundary layer by the two streamlines. The change 
in entropy at the edge of the boundary layer is determined from the entropy distribution 
along the shock. When the entropy at the edge of the boundary layer and the static pres- 
sure distribution over the body are known, the remaining edge conditions and flow proper- 
ties can be found from equilibrium air thermodynamic tables. The variation in boundary- 
layer edge conditions due to change in entropy may be calculated by balancing the mass 
flow in the boundary layer with that entering the stream tube through the shock given by 



(18) 


The expression for the shock radius is then given by 
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The value of the entropy at the edge of the boundary layer is found with rgy^r^ (from 
eq. (19)) from a table of shock radius and shock entropy from the inviscid calculation. 


Skin- Friction Calculation 


A value of the local coefficient of skin friction is needed at each point of calculation 
in order to evaluate the right-hand side of equation (3). In addition to its use in the solu- 
tion of the integral momentum equation, the local skin-friction coefficient may be used to 
obtain local heat transfer by using the Reynolds analogy to obtain a local Stanton number. 
In reference 30 a comparison has been made of four turbulent sMn-friction theories with 
experimental data. The results of the comparison indicated that the skin-friction laws 
of Van Driest II (that is, Karman mixing length, ref. 31) and Spalding-Chi (ref. 7) gave 
reasonably good agreement with heat-transfer data over 0.1 < T^/T^ ^ 1.0 when the 

( N«t e\ 

Reynolds analogy factor ^that is, j 1-0 for the Van Driest II theory and 1.16 


for the Spalding-Chi theory. The calculation of the Van Driest n theory in reference 30 
(described in detail in ref. 6) uses a Karman-Schoenherr incompressible skin-friction 
expression. The Karman-Schoenherr equation uses a transformed value of the local 
Reynolds number based on momentum thickness to obtain skin-friction coefficient for 
the incompressible plane. The skin-friction coefficient for the incompressible plane is 
transformed to the compressible plane by using an expression which is a function of Me, 
and T-w/Te, and is analogous to the Spalding-Chi Fq function. (See ref. 7.) This anal- 
ogous Fg function is referred to as the "Van Driest n Fc function" and is given in 
reference 6 as 


Fc = 


0.176Me2 
sin“^Q! + sin"^/3 


( 20 ) 


where a and /3 are functions of Mach number and temperature ratio. (See ref. 6.) 

It should be noted that the Van Driest n and the Spalding-Chi skin-friction theories 
presented in reference 6 use the same expression to obtain the incompressible skin fric- 
tion and also use the same F^ function (eq. (20)) to transform the incompressible skin 
friction to the compressible plane. The F^ function given in reference 6 uses the lin- 
ear Crocco in its derivation ^that is, Aj^ = 0, Br = 1.0, = 1.0 in eq. (10)^. The 

skin-friction calculations maide in this report by using the Spalding-Chi theory used the 
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Fc function as given in reference 7. The general expression given in reference 7 for 
Fc is 



The F(. function given in equation (21) may be integrated by using a real- gas vari- 
ation of equation (12); however, the correct values of otr must be used. The ideal gas 
F^ function used in reference 7 to correlate the skin-friction data is given by 



where Nj., as used in reference 7, is equal to 0.89. Equation (22) was derived from 
equation (21) with the assumption of a linear Crocco relationship (that is, = l.O). 

The local skin-friction coefficient obtained from the turbulent-flat-plate theories from 
either reference 7 or reference 6 must be corrected for application to a cone by the 
Mangier transformation factor (See appendix.) The local skin friction on a cone 

downstream of the end of transition is found from 

(t) = ^mt(t) 

V 2 /Cone \ ^ 'Flat plate 

where for fully turbulent flow 

Fmt “ 1-176 (24) 

The determination of the local value of skin friction from the Spalding- Chi (ref. 7) 
and Van Driest n (ref. 6) theories is calculated by use of a Reynolds number based on 
momentum thickness rather than on a Reynolds number based on a surface distance. The 
use of Re^e to find Cf/2 (1) eliminates the necessity of establishing a virtual origin 
of the turbulent flow and (2) makes use of the momentum thickness which is calculated 
directly from the momentum integral equation (eq. (3)). It is recognized that a pressure 
and entropy gradient can affect C^; however, this effect has been neglected in the pres- 
ent work. 
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Transition Region 


The value of Cf/2 in the transition region is assumed to vary as 

+ Y \ (Ci] _ /Cf\ ' 

2 2 2 tanh xj/ \ 2/turb \ ^/tr 

where 

Y = tanh ;// + tanh X 


( 25 ) 


X = xl/ - 2\l/ 

and \p = Constant. The value of xj/ controls the degree of nonlinearity of Cf/2 with 
respect to x/r^ in the transition region and for all calculations the value xp was equal 
to 2.0. The variation of Cf/2 in the transition region as calculated from equation (25) 
resulted in S- shaped curves which are characteristic of variations in body Cj/2 and 
heat- transfer data in transition regions. The initial value of skin friction at the start of 
transition may be taken from a solution for laminar flow or from a value based 

on experimental data. The value of skin friction at the end of transition (^f/^)j-^j.j 3 
calculated from the turbulent value found from equation (23). 

The value of N in the transition region is assumed to vary as 



N = N,r + rt^Nrb - Ntr) (26) 

where Y is the expression used in equation (25) and xp = 2.0. The value of N^j. is 
taken from the theoretical laminar calculation. The value of N at the end of the transi- 
tion Nturb is calculated from equation (8) . 

The value of ajji used in equation (12) varies linearly from an initial value at the 
start of transition to ajjj =1.0 at the end of the transition region. 


Heat Transfer 

The local Stanton number is calculated from a modified Reynolds analogy in the form 

i^St,e = ^ ^AF (27) 

where is the Reynolds analogy factor. The value of the Raf depends on which 

local turbulent skin-friction theory is used. If the Van Driest n (ref. 6) is used for heat 


13 



transfer, the Reynolds analogy factor is a function of h^He as listed in table I and for 
the Spalding-Chi skin-friction theory (ref. 7), Rat? = ^ . 


TABLE I.- VARIATION OF REYNOLDS ANALOGY FACTOR WITH hw/Hg 


hw/He 

Raf 

“^<0.2 

1.0 



0.2 ^ 0.65 

0.8311 + 0.9675 ^ - 0.6142&f 

He 

Hg \Hgy 

^>0.65 

1.2 

He 



The variation of the Reynolds analogy in table I is based on the results of refer- 
ences 32, 33, and 30 which indicated a decrease in R^f ^ decrease in hw/He- 

The results in reference 32 show considerable scatter in the data; however, there is a 
definite decrease in the value of the Raf with a decrease in h-^/He based on a least- 
squares curve fit of the data. The value of Raf the experimental data of refer- 
ence 33 showed a marked decrease for a ratio below 0.3. In reference 30 

Raf ® 1-0 (Nr = 0-9) was found for a simultaneous measurement of skin friction and heat 
transfer. Additional heat- transfer data in reference 30 indicate an increase in Raf 
as Tw/Tt increases similar to those of references 32 and 33. When the Spalding-Chi 
(ref. 7) theory was used in equation (27), Raf ~ (that is, Colburn's modified 

form of Reynolds analogy) was found to give the best prediction of heat transfer. The 
value of Npr is taken from tables of Prandtl number given in reference 29 evaluated at 
h' (see eq. (15)) and the local static pressure. 

The local heat-transfer coefficient h is calculated from 

h = Ng^^gUgP0 (28) 

The heating rate q is calculated from 

q = h(haw " M (29) 

where haw = he + Nr(Hg - hg) where Nr = 0.89. 

A constant value of Nr = 0.89 through the transition region is not correct; how- 
ever, for cold wall calculations (that is, h^ negligibly small compared with He), the 
maximum error at the start of transition is approximately 4 percent. When the wall 
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temperature becomes significant compared with the adiabatic wall temperature, the vari- 
ation of Nr should be included in equation (29). A summary of several of the parameters 
used in calculation of the transitional and turbulent boundary is given in table II. 

TABLE n.- SUMMARY OF PARAMETERS USED IN TWO REGIONS OF CALCULATION 


Parameter 

Initial value at 
start of transition 

Transitional boundary layer 

Turbulent bound^y layer 

N 

Laminar profile 

Equation (26) 

Equation (8) 

“m 

Laminar profile 

Linear variation from laminar 
value at start to 1 .0 at the 
end of transition region 

“m = 10 

Cf/2 

Laminar value (initial value 
is arbitrary; any value 
may be used) 

Equation (25) 

Equation (23) 

(Van Driest II (ref. 6) 
or Spalding- Chi (ref. 7)) 

Mangier transformation 
factor (Fjvix) 

Not applicable 

Not used (see eq. (25)) 

(Used in eq. (23)) 

Van Driest n — 1.176 
(eq. (33)) 

Spalding- Chi (eq. (34)) 

Reynolds analogy 
factor (Raf) 

The basis for the Raf same for the transitional and turbulent region, that is, 

Van Driest: Table I 

Spalding- Chi: 


METHOD OF CALCULATION 


The general method of calculation using the Langley Research Center computer 
program D-3340 is an iterative procedure which requires a calculation from the begin- 
ning of transition to the end of the body to be repeated until the velocity at the edge of 
the boundary layer changes less than 1 percent from one iteration to the next. When the 
velocity change is less than 1 percent, the calculation is terminated. The calculation pro- 
cedure first requires an inviscid solution (see refs. 27 and 34) followed by a viscous solu- 
tion which uses the wall conditions at the y = 0 boundary conditions and the results from 
the inviscid solution to establish the initial boundary conditions at the outer edge of the 
boundary layer. The conditions at the outer edge of the boundary layer change as the 
entropy at the edge of the boundary layer decreases from the stagnation value as a result 
of the mass flow into the boundary layer. If sufficient mass enters the boundary layer, 
the entropy at the edge will approach the sharp- cone value. The radius of the stream 
tube at the shock is determined from mass flow considerations. (See eq. (19).) The 
entropy at the edge of the boundary layer is determined from a table of shock radius 
against entropy as calculated from the inviscid solution. The viscous solution is then 
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inserted between the body and the inviscid solution, the assumption being that the dis- 
placement thickness correction does not change the effective body shape. Since there is 
no computed change in the effective body shape, the original pressure distribution and 
shock shape from the inviscid calculation do not change as a result of the viscous calcu- 
lation. The flow properties at the edge of the boimdary layer are determined from real- 
gas thermodynamic tables by using the pressure distribution and entropy at the boundary- 
layer edge to enter the table. 

Thermodynamic Properties for a Real Gas 

The flow properties in the viscous and inviscid flow fields, such as Pg, Pg, Tg, 
Se, he, Ue, etc., are calculated from the equilibrium air thermodynamic computer sub- 
routine described by Lomax and Inouye for an inviscid solution in reference 27. During 
the calculation, thermodynamic data at a point are found by entering the real-gas subrou- 
tine with pressure and entropy. For calculations of density profiles through the boundary 
layer, the real-gas calculation is modified to allow the thermodynamic data to be found 
for a given enthalpy (or temperature) and pressure. The procedure in this case is to 
enter the real-gas subroutine with various estimated values of entropy and the local pres- 
sure until the value of entropy is found that yields the desired value of enthalpy (or 
temperature) . 


Calculating Procedure 

Prior to the boundary- layer calculation, the inviscid flow field is determined by the 
Lomax and Inouye blunt-body and method-of-characteristics programs. (See refs. 27 and 
34.) The inviscid solution gives the first-order stagnation entropy flow property distri- 
bution along the body which is used as the initial conditions at the edge of the boundary 
layer. In addition, the shock shape rg/Tn entropy distribution along the shock are 
found from the inviscid solution. The results of the inviscid solution are first used to 
make a laminar boundary- layer calculation, with variable entropy, over the entire length 
of the body. (See ref. 35.) The turbulent boundary-layer calculations then follow. As 
previously noted, the pressure distribution over the body does not change as a result of 
the viscous calculation; however, all other conditions at the edge of the boundary layer 
change because of the change in local entropy. The turbulent-boundary-layer calculation 
is initiated at the point where transition has been determined to start. The initial edge 
conditions used in the turbulent-boundary-layer calculation come from the laminar solu- 
tion. The initial use of the edge conditions from the variable -entropy solution for laminar 
flow enables the turbulent calculation to be completed in less time than the initial use of 
the stagnation entropy edge conditions directly from the inviscid solution. Specifically, 
the output from the solution with variable entropy for laminar flow gives tables of rs/^n> 


16 


and dimensionless shear as a function of x/rn which are used for input to the turbulent- 
boundary-layer calculation. An initial entropy distribution at the edge of the boundary 
layer used in the first iteration for turbulent flow is determined from the table of rg/rn 
against x/r^ from the solution for laminar flow and the inviscid entropy distribution 
along the shock. At the point where the laminar-boundary- layer calculation stops and 
the transitional- boundary- layer calculation starts, the dimensionless shear table is used 
to find the initial value of Cf/2 and the rg/r^ table is used to find an initial value of 
rgy'rn. The value of Cg/rj^ at the start of transition specifies a value of entropy; with 
this value of entropy and the local static pressure, the flow conditions at the edge of the 
boimdary layer are known by entering the real-gas subroutine with entropy and pressure. 
The initial values of the velocity and density profile are based on the laminar profiles 
from which values of N and otjjj are determined, just upstream of the start of transi- 
tion. Nominally, the initial value of N is 1.5 and o!jjj is 1.0. Downstream of the end 
of transition, o!jjj is held constant (usually a value of 1 .0) . The value of the boundary- 
layer thickness at the start of transition is calculated from equation (19) by use of the 
initial value of rg^rj^, N, and Qijjj- The initial values of d/r^^ and S^J/r^ are found 
from equations (5) and (4), respectively. Downstream of the initial point of calculation, 
the value of S/r^ is calculated from the integration of equation (3), and 6/rn and 
sy'rj^ are calculated from equations (5) and (4), respectively. 

Transition Region 

The region from the end of the laminar -flow solution (beginning of turbulent solu- 
tion) to the point where the solution is considered to be fully turbulent is defined as the 
transition region. A guide for the selection of the length of the transition region, based 
on experimental data, is that the transition region ends when the local Reynolds number 
is approximately twice its value at the beginning of the transition region. (See ref. 36.) 
However, the guideline of a factor of two in Re,x transition is not always 

applicable when the edge conditions are changing very rapidly in the transition region 
because of the variable -entropy effects. These effects cause an increase in Reynolds 
numbers in the transition region which may be considerably different from a factor of 
2.0. Therefore, the location of the end of transition was put into the calculation as a 
specific value. The calculation procedure in the transition region is the same as the 
fully turbulent solution with three exceptions (see table n) : (1) the value of N in the 
transition region is calculated from equation (26); (2) the value of Cf/2 in the transi- 
tion region is determined from equation (25); and (3) the value of o!jjj varies from an 
initial value at the start of transition to value of 1.0 at the end of transition. After each 
calculation (iteration), from the start of transition to the end of the body, the relationship 
in the transition region for N and C^2 as a function of x/Tj^ (eqs. (25) and (26)) is 
recalculated for use in the next iteration. 
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Computer Calculation Time 


The time required for a given variable-entropy calculation for turbulent flow 
depends on (1) the number of iterative calculations from the start of transition to the end 
of the body; (2) the nvtmber of points on the body per iteration; (3) the number of points in 
the integration from the wall to the edge of the boundary layer; and (4) the error criterion 
for various iterative solutions used in the calculation. It usually took four iterations to 
find a reasonable engineering solution which had a change in velocity of less than 1 per- 
cent < O.Olj from one iteration to the next. For the flight calculation, the Runge- 

Kutta integration was limited to maximum step size of 16.0; however, if calculation 

^n 

time on the computer is extremely critical, the step size could possibly be increased 
with little effect on the final solution. The number of points required in the Gaussian 
quadrature for the calculation of S/fn 6*/rn from equations (5) and (4), respec- 
tively, was determined from a comparison of the values of 6/rn and calculated 

for various numbers of points used in the quadrature. It was determined that a minimum 
number of 20 points could be used in the integration through the boundary layer. The 
momentum integral equation (eq. (3)) is solved by a variable -step- size fifth-order Runge- 
Kutta numerical scheme which uses a maximum relative error criterion of 0.001 for the 
value of 9 for one step of integration. The maximum percent error allowed in the itera- 
tive solution for the Spalding- Chi (ref. 7) Cf/2 was 0.4. When an iterative procedure 
had to be used in the determination of thermodynamic properties from the real- gas sub- 
routine (see ref. 27), the maximum allowable relative error was 0.001. For real-gas 
turbulent-boundary- layer calculations made at altitudes from 18.29 km to 25.91 km 
(60 000 ft to 85 000 ft) at M « 20, the time per point of calculation was approximately 
2.6 seconds on the Control Data 6600 computer system based on a single iteration from 
the start of transition to the end of the body and using the relative error criteria. 

RESULTS AND DISCUSSION 

Comparison of Ground Test Data and Theoretical 
Predictions of Heat Transfer 

In order to evaluate the integral theory in terms of its ability to predict heat trans- 
fer, a comparison was made between theoretical and experimental results for three sets 
of heat-transfer data. 

Flat-plate da ta.- In figure 2 six integral calculations of Stanton number were made 
in order to evaluate the effect of 0 !j^ and F(. (see eq. (12)) on the integral method when 
the skin-friction laws of Spalding-Chi (ref. 7) and Van Driest n (ref. 6) are used. Four 
calculations using the Spalding-Chi skin-friction theory were made for two Fq functions 
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X, cm 

Figure 2.- A comparison of heat-transfer data on a flat plate with integral and 

T 

finite-difference methods of prediction at Mach 6.0 and — = 0.2. 


(eqs. (21) and (22)) with equal to 1.0 and 2.0 for both methods of obtaining F(>. The 

other two calculations were made with the Van Driest n skin-friction theory, by using 

equation (20) for the Fq function, with values of 0 !p{ of 1.0 and 2.0. The six integral 

T 

calculations are compared in figure 2 with the heat-transfer data at of refer- 

ence 37. Five of the theoretical predictions, which use a linear Crocco relationship to 
derive the Fq equation, agree well with the data at peak heating but tend to overpredict 
by approximately 10 percent at the 50.8-cm (20-in.) station. The sixth Stanton number 
calculation from the Spalding- Chi theory, which uses ap = 2.0 to integrate the basic 
expression for Fg (eq. (21)), predicts a Stanton number which is approximately 
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70 percent above the data. In general, the results of the integral calculations show that 
when the function is derived for a linear Crocco relationship (that is, ap = 1.0, 
eqs. (20) and (22)), the heat-transfer prediction is close to the data even though the inte- 
gral parameters are calculated for 0 !j^ = 2.0. Therefore, if equation (21) is to be inte- 
grated for either an ideal gas or a real-gas variation of density through the boundary 
layer, evidently an 0 !ju = 1.0 must be used. In addition, the predictions calculated from 
the finite-difference theory of reference 38 for both the laminar and turbulent flow are in 
general agreement with the five integral predictions. 

Blunt-cone data.- Two sets of heat-transfer data are used as a basis of comparison 
for heat-transfer predictions with variable entropy and are shown in figures 3 and 4. The 



Pigtire 5-- A comparison of normalized heating-rate predictions from the final 
iteration's for laminar and turbulent boundary layers to experimental wind- 
tunnel data at = 11.5 on a 5° half-^gle cone at R = 8.5 X 106/m 
(2.6 X 10^/ft). r^ = 2.54 imn (O.IO in.). 
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data from reference 39, shown in figure 3, were measured in a "shock tunnel" at a free- 
stream Mach number of 11.5 on a 3.66-meter-long (12-ft), 5° half-angle cone, with a 
0.254-cm (0.10-in.) nose radius. The data (ref. 39) in figure 3 were normalized in 
order that the laminar data match the theoretical prediction for laminar flow (ref. 35) at 
X = 50.8 cm (20 in.). When the data from reference 39 are normalized, the theoretical 
prediction in figure 3 for both laminar flow (ref. 35) and the integral theory for turbulent 
flow, presented in this paper, agree reasonably well. 

A comparison between data obtained by P. Calvin Stainback of Langley Research 
Center and theoretical predictions from two integral methods, using Spalding-Chi (ref. 7) 
and Van Driest II (ref. 6) are shown in figure 4. The data were taken on a 73.66-cm-long 
(29-in.), 5® half-angle cone with 0.508-mm (0.020-in.) nose radius. The tests were at a 
free- stream Mach number of 8.0 and a unit Reynolds number of 5.44 x 107 per meter 
(1.658 X 107 per ft). The heating-rate data in the region of turbulent flow in figure 4 is 
considered to be reliable because of the accurate measurements of total temperature, 
wall temperature, and skin thickness, and also because of the proven data reduction 
method for obtaining the heating rate from the rate of change of the wall temperature. 

(See ref. 40.) The cone used in the tests had a circumferential rearward-facing step 
0.254-mm (0.010- in.) high at a station 3.81 cm (1.5 in.) from the nose. The purpose of 
the step was to simulate a flight configuration. The combination of the rearward-facing 
step and the uncertainty in the skin thickness in the small nose piece are believed to be 
the reasons that the laminar theory (ref. 35) slightly underpredicts the data upstream of 
the 15.24-cm (6.0-in.) station. In figure 4 the calculation of heat transfer in the region 
of peak heating for fully turbulent flow which uses the Spalding-Chi skin-friction theory 
(ref. 7) slightly overpredicts the data; whereas in the same region the Van Driest n (ref. 6) 
theory slightly underpredicted the data. The difference in the prediction in peak heating 
is a result of the two theories using different variations of Fjyf-p in the turbulent flow 
region. The Van Driest n skin-friction theory used Fjy^fp = 1,176 over the entire turbu- 
lent flow region whereas the Spalding-Chi theory uses equation (34) (see appendix) to 
obtain the variation of Fj^>p. At the end of the cone, both skin-friction laws use 
Fmt = 1-176 and both theories are in reasonably close agreement. The R^f ’^sed 
with the Van Driest II skin-friction theory is given in table I, and the Spalding-Chi skin- 

rt /o 

friction theory uses R^p = (Npj.j ' . 

Variable Entropy Turbulent Boundary Layer - Sample Flight Calculations 

The results of a real-gas variable-entropy turbulent-boundary- layer calculation 
made for sample flight conditions are shown in figures 5 to 17. The sample calculations 
herein were made for a blunt 5° half-angle cone with a 1.016-cm (0.4-in.) nose radius 
at an altitude of 18.29 km (60 000 ft) at a velocity of 5.80 km/sec (19 019 ft/sec). The 
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Figure 5 -- The effect of three entropy distributions at the edge of the bovindary 
layer on the local Mach number for three turbulent -boundary- layer calcula- 
tions at 18.29 hm (60 000 ft) and Uoo = 5-80 km/sec (19 019 ft/sec), 
r^ = 1.016 cm (0.4 in.). 

sample calculations for the turbulent boundary layer on a 3.96-meter-long (13-ft) cone 
were made with the assumption that the boundary layer was nearly turbulent over the 
entire body. However, in order to simulate the actual method of calculation, a short dis- 
tance of laminar boundary layer was assumed to precede the start of the turbulent calcu- 
lation. For the illustrative sample calculation, the start and end of transition were 
assumed to be 15.0 and ^=^3^, respectively. The assumption of the extremely 

long transition region was based on the combined effects of (1) bluntness and (2) the varia- 
ble entropy. Both of these effects are accentuated by the far upstream location of the 
start of transition. The turbulent-boundary-layer calculations in figures 5 to 15 show 
the results for three sets of flow properties at the edge of the boundary layer, hereafter 
referred to as edge conditions. The three edge conditions used in the calculation resulted 
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Figiare 7.- The effect of three entropy dlstrihutions at the edge of the boundary 
layer on (l) the ratio of the local static enthalpy to the inviscid flow field 
total enthalpy and (2) the ratio of local velocity to free-stream velocity for 
three turbulent -boundary- layer calculations at 18.29 km (60 000 ft) and 
i 3 co = 5-80 km/sec (19 019 ft/sec), r^ = I.016 cm ( 0.4 in.). 


from the following three entropy distributions at the edge of the boundary layer: (1) con- 
stant stagnation entropy; (2) the entropy from the final laminar solution (see ref. 35); and 
(3) the entropy resulting from the final turbulent calculation. For the variable -entropy 
calculations for a turbulent boundary layer described in this paper, the calculation for 
the first iteration uses the entropy distribution resulting from the final iteration for a 
variable-entropy solution for laminar flow instead of using the stagnation- entropy distri- 
bution. The stagnation- entropy results are presented to illustrate the results that would 
be obtained if variable entropy were not considered in the calculation. 
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Figure 8.- 'The effect of three entropy distributions at the edge of the 
boundary layer on the value of the IT po-wer-law exponent for three 
turbulent-boundary- layer calculations at 18.29 hm (6o 000 ft) and 
Uj^ = 5.80 km/sec (19 019 ft/sec), rn = I .016 cm (0.4 in.). 

Flow properties at edge of boundary layer .- The effects of the three entropy distri- 
butions on the Mach number, the Reynolds number, and the h^Hg ratio at the edge of 
the boundary layer are shown in figures 5,6, and 7, respectively. The variation of the 
Mach number at the edge of the boundary layer, shown in figure 5, ranges from a value 
of approximately 3.6 at the upstream end of the cone to a downstream value of approxi- 
mately 15.1 for the final iteration for turbulent flow. Thus, at these flight conditions a 
sharn-cone value of Mp = 14.6 is reached in the final calc ulation for turbulent flow at 
an x/rn of approximately but downstream of ^ the Mach number overshoots 

the sharp- cone value of 14.6. The overshoot in the local Mach number is a result of low- 
entropy air (that is, less than the sharp-cone value) entering the boundary layer in a region 
where the surface pressure is at approximately a sharp- cone value. If the cone were long 
enough, the local value of Mach number would decrease gradually and become asymptotic 
to the sharp-cone value of 14.6. 

The effect of variable entropy at the edge of the boundary layer on the local Reynolds 
number based on x/r^ distance is shown in figure 6. The value of Rg^x from the final 
calculation for turbulent flow at the end of the cone increases by a factor of approximately 
72 compared with the value found for the stagnation entropy at the same location. 
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Figure 9 *- The effect of three entropy distrihutions at the edge of the houndary 
layer on the momentum thickness for three turbulent -boundary-layer calcula- 
tions at 18.29 km (60 000 ft) and u^^ = 5*80 km/sec (19 019 ft/sec), 
r^ = 1.016 cm (0.4 in.). 

JlpS 

The Re x the final turbulent calculation is equal to the sharp-cone value at — ~ ■8 66 - 
’ liflS Tn 

and for ^ > - 36 5 ‘ the value of Re,x overshoots the sharp- cone value as was found in fig- 
ure 5. Again if the cone were long enough, the value of Re x from the final calculation 
would approach the sharp- cone value. 

The effect of variable entropy on the local static enthalpy and local velocity is shown 
in figure 7. The he/He ratio at the end of the cone for the final turbulent calculation is 
approximately 20 times less than the h^He ratio calculated from the stagnation- entropy 
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Figure 10.- The effect of three entropy distributions at the edge of the boundary 
layer on boundary- layer thickness for three turbulent -boimdary-layer calcula- 
tions at 18.29 hm (60 000 ft) and = 5 . 8 O km/sec (19 019 ft/sec), 
rji = 1.016 cm ( 0.4 in.). 


edge condition. The Ue/u^o ratio at the end of the cone for the final turbulent calculation 
is approximately 1.35 times greater than the Ug/u^o ratio calculated from stagnation- 
entropy edge conditions. The hg/Hg and Ug/u^ ratios from the final calculation for 
the^^rbulent boundary layer approach the sharp-cone values at approximately an x/r^ 
of 80^ The values of hg/Hg, Pg, Tg, Mg, and Ug are shown in table m, at 
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TABLE III,- THE VALUES OF VELOCITY, DENSITY, TEMPERATURE, 
AND MACH NUMBER AT #- = 400 FOR STAGNATION 
AND VARIABLE ENTROPY CONDITIONS 
[Altitude, 18.29 km (60 000 ft); = 5.80 km/sec (19 019 ft/sec^ 



Stagnation 

entropy 

Final turbulent 
calculation 

he/He 

0.4655 

0.02133 

Pg, kg/m3 (Ibm/cu ft) 

0.0314 (0.00195) 

0.4019 (0.02509) 

Tg, K (Or) 

3960 (7128) 

358.7 (645.6) 

Mg 

3.34 

15.15 

Ug, km/sec (ft/sec) 

4.264 (13 990) 

5.77 (18 930) 


= 400, for an entropy distribution change from a stagnation value to the final turbulent 
value. The flow properties of h^Hg, Pg, Tg, Mg, and Ug vary by factors of 21.8, 
12.9, 11.0, 4.5, and 1.4 respectively. Thus, for the same change in entropy distribution, 
the variation of flow properties can be as small as a factor of 1.4 (ug) and as large as a 
factor of 21.8 (h^Hg). 

Bound ary- layer parameters .- The boundary- layer parameters ^N, 0/rj^, 6/'rj^, 

and 6^rn) shown in figures 8 to 11 were calculated from the same boundary -layer cal- 
culation used to determine the edge properties shown in figures 5 to 7. The variation of 
the N "power-law” exponent for the three sets of conditions at the edge of the boundary 
layer is shown in figure 8. The value of N at the end of the cone = 400j varies 
from approximately 4.3 for the edge conditions based on stagnation entropy to a value of 
approximately 9.8 for the final calculation for a turbulent boundary layer. The calcula- 
tions for all three entropy distributions start at the beginning of transition with N = 1.5 
based on a laminar boundary- layer profile. 

The dimensionless momentum thickness shown in figure 9 for the three edge 

condition distributions is calculated from an integration of equation (3). The decrease in 

0/rj, from the start of the final iteration to an x/rj, of approximately 100 would imply 

d0/rn 

a negative skin friction, that is, for sharp-cone conditions where — — — is equal to 

dx/rn 

Cf/2. However, the reason for a decrease in is not from a negative value in Cf/2 

but from the combined effects of variable entropy and transitional flow which is felt in 
equation (3) in terms of: (1) the value of the edge condition gradients; (2) the variable- 
entropy term (eq. (6)); (3) the value of term; and (4) the value of Cf/2, 

Thus, the combination of these effects has a significant effect on the integration of the 
momentum equation and the resulting value of 6/t^. 
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PigTxre 11.- The effect of three entropy distributions at the edge of the boundaiy 
layer on the displacement thickness for three turbulent -boundary- layer calcu- 
lations at 18.29 km (60 000 ft) and u„ = 5.8O km/sec (19 019 ft/sec), 
rn = 1.016 cm (O.if in.). 

The effect of three different edge conditions on the calculation of the boundary- layer 
thickness 6/rn the displacement thickness S^yrn is shown in figures 10 and 11. 
These results show no common trend in the growth of 5 / tj ^ and 6 with respect 
to the variation of the three edge (entropy) conditions. However, it appears because of 
the rapidly changing edge conditions, resulting from a change in local entropy from a 
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Figure 12.- The shock radius distrihution based on a mass rate of flov into the houadaiy 
layer for two turbulent- and one larainar-boundary-layer calculations at 18.29 km 
(6o 000 ft) and Uco = 9-80 km/sec (l9 019 ft/sec), r^ = 1.0l6 cm (0.4 in.). 


stagnation to a sharp-cone value, that and 5*/rj| may not follow the classical 

variation with ( that is, ^ a 1-_ and ^cc but are strongly Mach 

' R 1/5 


R* 


1 / 5 , 


V ,x ,x 

number dependent, the value of and 6*/rn increasing with increasing Mach 

number. 


The shock radius Vgy'rj^ (see eq. (19)) shown in figure 12 is used to determine the 
value of entropy at the edge of the boundary layer. The values of rg/r^ in figure 12 
labeled "Final laminar iteration" were calculated during the final iterative solution for 
laminar flow (ref, 33) and were used to determine the entropy distribution for the first 
iteration in the solution for turbulent flow. The rg/rn values calculated for the turbu- 
lent bovmdary layer from stagnation entropy distributions are greater than the values cal- 
culated from the final iteration for a laminar boundary layer; this result indicates that 
even for a low Mach number (that is. Mg « 3.5), the turbulent boundary layer has a larger 
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mass entrainment than the laminar boundary layer. It can be seen from figure 12 that the 
mass entrainment for the final calculation of the turbulent boundary layer is very large 
compared with the calculation of mass entrainment for either the laminar boundary layer 
or the turbulent boundary layer with stagnation conditions at the edge of the boundary 
layer. 

Skin friction .- The local skin friction C^2 used in the calculations for figures 5 
to 12 is shown in figure 13, The skin friction is calculated from the Van Driest n (ref. 39) 
theory, (See table n,) The initial value of Cf/2 at the start of transition is taken from 

the solution for laminar flow ^that is, ^ = 1.02 x 10"3j. For the two "first iteration" 
calculations, equation (25) cannot be used and, as a result, the value of Cf/2 jumps from 
the initial laminar value to the value calculated by the Van Driest n (ref. 6) theory. The 
level of Cf/2 in figure 13 decreases as the entropy distribution changes from (1) stag- 
nation, (2) to laminar calculation, and (3) to a final iteration value. In terms of the clas- 
sical variation of Cf/2 with Re^xj decrease in Cf/2 for the final iteration calcu- 
lation would be expected, based on the corresponding rapid increase in the value of Rq x 
shown in figure 6. However, the variation of Cf/2 for the final calculation is more com- 
plex, in terms of shear and heat transfer, than would appear from the results in figure 13. 

The variation of the wall shear calculated from 


'w 


Cf 


(30) 


is shown in figure 14 for three entropy distributions. The shear for the final calculation 
for turbulent flow shows a rapid rise at the start of transition to a nearly constant value 
for the fully turbulent flow. However, the rapid rise in the transition region is not 
reflected in the corresponding value of Cf/2 in figure 13, The reason for the rapid rise 
in shear can be seen by examining the variation of PgUg2 (fig. 15) through the transition 
region as calculated for the same edge conditions used in figures 13 and 14, A compari- 
son of the final iteration results in figures 14 and 15 from the beginning to the end of tran- 
sition region shows that PgUg2 increases by a factor of approximately 25 whereas 
only increases by a factor of 5. These increases result in a factor of 5 decrease in Cf/2 
in the transition region. It should be noted that the heat transfer in the transition region 
is similarly affected by a rapid rise in the product of p^Ug, The rise in p u when 
coupled with a nearly constant value in h^^ - produces a rapid rise in the heating 
rates in the transition region. (See eqs, (28) and (29),) 

Effect of nose radius .- The effect of nose radius on the local Mach number and 
Reynolds number is shown in figures 16 and 17, respectively. The results shown in these 
two figures were taken from the final iteration of a variable-entropy calculation for tur- 
bulent flow at zero angle of attack, at an altitude of 18.29 km (60 000 ft), and a velocity of 
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Figure 13 •- The effect of three entropy distributions at the edge of the boundary 
layer on the local skin- friction coefficient for three turbulent -boundary-layer 
calculations at 18.29 km (60 000 ft) and u^o = 5-80 km/sec (19 OI9 ft/sec). 

= 1.016 cm (0.4 in.). 
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End of transition 



Final iteration 


^ First iteration 

(final laminar entropy distribution) 


First iteration (stagnation entropy) 


Figure l 4 .- The effect of three entropy distributions at the edge of the boundary 
layer on the shear stress at the -wall for three turbulent -boundary-layer cal- 
culations at 18.29 km (60 000 ft) and = 5-80 km/sec (19 019 ft/sec), 
rn = 1.016 cm ( 0.4 in.). 
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End of transition 



Figure 15 •- The effect of three entropy distributions at the edge of the boundary 
layer on the local value of Pe'^e^ three turbulent -boundaiy- layer cal- 

culations at 18.29 km (60 000 ft) and u^o == 5*80 km/sec (19 019 ft/sec), 
rn = 1.016 cm (O.^l- in.). 
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Figure l6.- The effect of nose bl-untness on the local Mach number for the final 
iteration calculation on the turbulent boundary layer at 18.29 kni (60 000 ft) 
and = 5.80 km/sec (19 019 ft/sec). 
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Figure 17 •- The effect of nose hluntness on the local Reynolds number for 
the final iteration calculations on the turbulent boundary layer at 
18.29 3nn (60 000 ft) and u__ =5-80 tan/sec (19 019 ft/sec). 



5.80 km/sec (19 019 ft/sec) for nose radii of 1.016 cm (0.4 in.) and 0.399 cm (0.157 in.). 
The start of transition was again assumed to be at x = 15.24 cm (6.0 in.), and the end of 
transition was assumed to be at 2.26 m (89.0 in.) for both nose radii. The results in fig- 
ures 16 and 17 show that the smaller nose radius (r^ = 0.399 cm (0.157 in.)] approaches 
a sharp- cone value of Mach number and Reynolds number faster than the large nose radius 
(rji = 1.016 cm (0.40 in.)]. For both nose radii the local Mach number and local Reynolds 
number overshoot the sharp-cone values. However, at the end of the cone that has a 
0.399 cm (0.157-in.) nose radius, the values of Re,x ^e h^ve become asymptotic 

to the sharp-cone values. 

CONCLUDING REMARKS 

A simple integral technique for calculating the transitional and turbulent bovindary 
layers for equilibrium air on an axisymmetric body at zero angle of attack - including 
the effects of variable entropy — has been presented. Boundary- layer velocity profiles 
were calculated from an N power-law correlation. Density profiles were calculated 
by use of a modified Crocco relationship. Skin friction was determined from either the 
Van Driest n or Spalding-Chi theories, by use of a momentum thickness Reynolds num- 
ber. Heat transfer was calculated by using a modified form of Reynolds analogy. The 
Reynolds analogy factor, used with the Van Driest II skin-friction theory, varied as a 
function of wall enthalpy ratio h^/Hg from 1.0 to 1.2, but with the Spalding-Chi theory 

it was held constant at (Npr)"^^^. 

Calculations of heat transfer for transitional and turbulent flow agreed with data 
measured on a flat plate at Mach 6.0 and two blunt cones at Mach 8.0 and 11.5. The cal- 
culations for the flat plate showed that the transformation function Fg used to transform 
the skin friction to the compressible plane must be calculated with the assumption of a 
linear Crocco relationship. The predictions of heat transfer on the flat plate obtained by 
use of the integral theory are in agreement with a finite-difference theory. Both the 
Van Driest II and Spalding-Chi skin-friction theories gave a reasonable prediction of heat 
transfer for the Mach 8.0 cone data provided the proper Reynolds analogy factor was used. 

Sample turbulent- boundary- layer flight calculations were made at an altitude of 
18.29 km (60 000 ft) and flight conditions of approximately 6.09 km/sec (20 000 ft/sec) 
for a 3.96-meter-long (13-ft), 5° half-angle cone with a nose radius of 1.016 cm (0.40 in.). 
The calculations were made for three entropy distributions; (1) stagnation value; 

(2) laminar-boundary-layer values; and (3) turbulent-boundary- layer values. The cal- 
culation showed that mass entrainment was considerably larger for the turbulent-boundary- 
layer distribution than for the other two distributions; as a result, a sharp-cone entropy 
was reached at approximately 410 nose radii from the nose of the cone. The calculation 
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showed that as the edge entropy decreased, the shear at the wall increased considerably, 
whereas at the same time the skin-friction coefficient showed a marked decrease. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., April 6, 1971. 
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APPENDIX 


TRANSFORMATION OF THE LOCAL AND AVERAGE TURBULENT SKIN 
FRICTION FROM A FLAT-PLATE TO A SHARP-CONE VALUE 


Local Skin Friction 

The general form of the Mangier transformation factor Fj^-p (see ref. 41) con- 
verting local flat-plate skin friction to a sharp-cone value is 


Fmt = 


Cf 

Cone _ 

JL 

rn 


n 

n-1 

Cf 

Flat plate 

Jo 1 




(31) 


When the Blasius form of the turbulent skin friction for a flat plate 

■ 1/n 


(32) 


Cf oc Re,x 

is used, with n = 5.0, along with the assumption of turbulent flow over the entire cone, 
the resulting Mangier transformation factor is 

Fmt = (33) 

In the region of fully turbulent flow, the Fjyj^ = (eq. (33)) was used for the 

calculations in which the Van Driest H skin-friction theory was used. When the Spalding- 
Chi skin-friction theory was used, an alternate variation of the F|^-p was used which 
took into account the transition from a laminar Fi^T to a turbulent F^t (1.176) 

in the form 


FmT^ “ 2V3Fmt + 3.0 = P 


X 

l^o] ' 

rn " 

Wq 


(34) 


where 


P = 


0.309136 


m 


max 


Xvo\ 


The ^XyQ^rn^^ term was evaluated at the beginning of fully turbulent flow from equa- 
tion (9) and the was evaluated at the end of the cone. 

A check was. made on the effect of nose bluntness on the value of Fj^t found from 
equation (31) for a 3.96-meter-long (13-ft), 5° half-angle cone with a 1.016-cm (0.4-in.) 


40 


APPENDIX - Concluded 


nose radius. The check showed that the value of which included the effects of 

nose bluntness was less than 1.0 percent smaller than the value of 1.176 (eq. (33)) found 
for a sharp cone. 


Average Skin Friction 

The ratio of the average skin friction on a sharp cone to that on a flat plate is 


Cone 


(Flat plate 




(35) 


where L is a constant and is the same for the flat plate and cone. When equations (32) 
and (33) are substituted into equation (35) and the indicated integration is performed, the 
resulting ratio of the average skin friction on a cone to that on a flat plate is 


Cxj! 

_ — .= 1.045 

Cp 

Flat plate 


(36) 
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